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We investigate the properties of a spinless Fermi gas close to a p-wave interaction resonance. We 
show that the effects of interaction near a p-wave resonance are captured by two contacts, which are 
related to the variation of energy with the p-wave scattering volume v and with the effective range 
R in two adiabatic theorems. Exact pressure and virial relations are derived. We show how the two 
contacts determine the leading and sub-leading asymptotic behavior of the momentum distribution 
(~ and ~ and how they can be measured experimentally by radio-frequency and photo¬ 

association spectroscopies. Finally, we evaluate the two contacts at high temperature with a virial 
expansion. 


Introduction. In the past decade, degenerate Fermi 
gases close to scattering resonances have attracted both 
theoretical and experimental attention [T] . In the unitary 
Fermi gas close to an s-wave resonance, it is understood 
that thermodynamic properties are universal [2] , depend¬ 
ing only on a single function, called the “contact” [H6]. 
Its manifestations in physical properties have been ex¬ 
tensively explored and confirmed in experiments 
Extension to arbitrary dimensions has been considered 
[Mill]. However, so far contact has only been consid¬ 
ered for s-waves, even though p-wave and higher partial 
wave resonances have been explored experimentally [H- 
[23] and theoretically [21H3Hj- 

In this Letter, motivated by the recent radio-frequency 
(rf) spectroscopic data near a p-wave Feshbach resonance 
in "‘'’K [39], we generalize the concept of contact to p- 
waves. In the case of an s-wave resonance, a single con¬ 
tact, which depends on the s-wave scattering length Us, 
is sufficient for the characterization of universal proper¬ 
ties of the system. For example, the two-body binding 
energy is given by where M is the mass of the 

atoms, and the effective range correction is in general 
small [33. In the case of p-wave scattering, however, the 
phase shift is given by cot 6{k) = —1/vk^ — 1/Rk for a 
short-range potential, and the effective range R is of fun¬ 
damental relevance, in addition to the scattering volume 
V. This can be seen clearly in the binding energy of a 
shallow p-wave bound state Et = h?R/Mv [31], depend¬ 
ing crucially on both v and R. As a result, to capture 
the universal properties of a spinless Fermi gas around a 
p-wave resonance, it is necessary to introduce two con¬ 
tacts, related to the variation of v and R, separately. We 
show how two adiabatic theorems [see Eqs. 0 and @] 
can be established and how the two contacts relate to the 
leading 1/^^) and the sub-leading (~ 1/fc^) terms of 
the high-momentum distribution. We also show how the 
two contacts can be measured spectroscopically. Finally, 
we use a virial expansion to determine each contact as a 
function of T, v, and R at high temperature. 


General Formulation. To start, let us consider the 
two-body problem, where two identical fermions of mass 
M interact via a short-range potential U{r) of range 
ro, tuned close to a p-wave resonance. The relative 
wave function in the p-wave channel can be written as 
^fc(r) = Xki'c)Yimif)/r, where m labels the projection of 
angular momentum along z-direction and k is the rela¬ 
tive wave vector. For low-energy p-wave scattering, the 
radial wave function Xk (c) can be expanded in powers of 
k^, Xk{r) = -|- -!-•••. In the asymptotic 

regime where 1/k r ^ tq, we fix the normalization 
such that the explicit form of Xfc(r) [and hence x^'^^ ^^nd 
X^^)] is 


Xk{r) 



+ .... ( 1 ) 


It is important to note that the above asymptotic forms 
for ^Iso hold for any shallow p-wave bound 

state in the corresponding asymptotic regime. Once the 
asymptotes are determined through Eq. Q, the short- 
range form (r < rg) of is completely fixed by 

two-body physics, due to competition between kinetic 
and potential energy and in particular, independent of 
the asymptotic wave vector k [5]. 

To proceed to the many-body case, we first need to 
derive two important identities, relating the change of 
V and R to that of the variation of the potential U{r). 
Consider two slightly different potentials U±{r) = U{r)± 
6U{r)/2, each with scattering volume v± and effective 
range R±. The radial Schrddinger equation is 


\ M dr"^ 


+ U±{r) + 



X±{r) 


h?k\ 

M 


X±{r). 


( 2 ) 


The term 2h?/{Mr"^) gives the p-wave centrifugal poten- 
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tial. Following the standard procedure Em], we find 
M 

6v~^ = --^J dr 6U{r)\x^°'>{r)\'^, (3) 

9M 

SR-'^ = — drSU{r)x^°\r)x^^\r), (4) 

where 3v~^ = — vZ^ and similarly for SR~^. 

Next, consider a spinless Fermi system of total 
number N confined in volume Zl with density n = 
N/H. = kp/{67r‘^) where kp is the Fermi wave vec¬ 
tor. The two-body density matrix p 2 (ri,r 2 ) = 

('0^(ri)'0^(r2)'0(r2)'0(i’i))) where '!/)l(r) creates a fermion 
at position r, is Hermitian and can be diagonalized 

P2{ri,r2) = ^naC(j'i,r2)(/'a(ri,r2). (5) 

a 

The eigenvalues ria satisfy the condition J2a = N{N— 
1), and the associated pair wave functions {<^a(ri,r 2 )} 
form an orthonormal set. In a rotationally invariant sys¬ 
tem, they can be further written as 


are identical to that of and x*'^^- Furthermore, 
for r < rg, both are uniquely fixed by the two-body 
physics. Thus, when evaluating the expectation value 
of any short-range function such as potential U{r), 
can be taken out of the integration over k. 

The interaction energy of the many-body system 
can be written in terms of p 2 as (U) = |/ C^(|ri — 
r 2 |)p 2 (ri,i" 2 )rf^rid^r 2 - Using the decomposition Eq. ([^ 
and Eq. 0 , we hnd 

^ E / drU\x^°'>\^ + J drUx^^h^^^], 

( 8 ) 

where we have defined two p-wave contacts for each 

m, 



\ J dk j dk'ajkajk'ik'^ + k'’^). (10) 

P,i 


'ii’a(ri,r2) = ^hexp(iP • R)(pji{r)Yem{r), (6) 
V Ur 

where R = (ri -|- r2)/2 is the center of mass and r = 
ri — r 2 is the relative coordinate and r = |r|. P can 
be regarded as the center-of-mass momentum of a pair 
and j, i, m label the quantum numbers of the relative 
radial direction, the angular momentum, and its z pro¬ 
jection respectively. Here the index a = {P,j,£,m} is 
a shorthand for all the quantum numbers that label the 
pair wave function. In a single-component Fermi gas, £ 
must be odd. In the region r —>• 0, Pji{r) ~ the 

p-wave channel has the strongest penetration inside the 
interaction potential U{r). As a result, we shall concen¬ 
trate only on the p-wave component, since it gives the 
dominant contribution to the interaction energy of the 
system. 

The pair wave function (/)Q,(ri, r 2 ), and hence ^Pji{r) is 
not an eigenfunction of the two-body Schrodinger equa¬ 
tion, but can be expanded in terms of the p-wave func¬ 
tions (setting £ = 1 and neglecting the subscript i from 
tpji thereafter) 

poo 

^i(^)= / dkajkXkir)aj^,Xt^{r), (7) 

Jo 

where {ajk-,ajK,} are the real expansion coefficients, the 
integration is over all scattering states, and we have also 
taken into account the possibility of a shallow bound 
state with radial wave function Xk('c) and binding energy 
Eb = Hk^/M = fi^R/{Mv), when v > 0 and R > 0. Ex¬ 
tension to multiple bound states is straightforward. An 
important consequence of such considerations is that, in 
the asymptotic region where rg r <C the form 
of Xfc(^)) s-nd hence Pjir), when expanded in power of 


Here the contribution from possible bound states is im¬ 
plicitly included in the integration over k. We note that 
has dimension of length, while has dimension 
of inverse length. Just as in the s-wave case, encap¬ 
sulate all the short-range correlations of the many-body 
system. As a byproduct, the two-body density matrix 
for r = ri — r 2 in the asymptotic regime rg ^ r ^ fcp ^ 
can be written as 



Now suppose that the potential U{r) can be con¬ 
trolled via an auxiliary parameter A, such that a small 
change in U{r) can be written as {dU/dX)d\. We can 
use the Hellmann-Feynman theorem and write dE/dX = 
{dZH/dX) = {dU/dX), where H = K,-\-U is the total many- 
body Hamiltonian with 1C denoting the kinetic energy, 
independent of A. Using Eqs. ( |3|4|8[ ), we find 


dE 


dv~^ 


R 


2M^ " 


dE 


’ dR-i 




( 12 ) 

In the simplest case of a shallow p-wave two-body 
bound state, the wave function in the asymptotic re¬ 
gion 1/k >> r :g> rg is given by '0K(r) = •\/A(l/r^ -|- 
K/r) exp{—Kr)Yim{r) with k = Rjv. It is then easy to 
obtain p 2 (r/ 2 , —r/2) = 2/H(i?/r'^ —ii^/ur^)|Yim(r)P for 
the two-body bound state m- One can extract directly 
that = 2R and = —2R‘^jv, consistent with 


the adiabatic theorems Eq. (12). The derivation also ap¬ 


plies to thermal equilibrium, in which case, one should 
replace the energy E by the free energy F of the system 
and keep temperature constant [^. 
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As in the s-wave case, a pressure relation and virial 
theorem can be found. In a uniform system, the uni¬ 
versal hypothesis is that the free energy can be written 
as F(T/Tf, kpV, kpR) close to a p-wave resonance, where 
Tp is the Fermi temperature. Using dimensional analy¬ 
sis [3H5], 


where i,j label the initial and final states, and pi de¬ 
notes the initial state distribution. In the region Ep <C 
hijj «C Ep = k? /MR^, one finds [421144] 


rrf(w) = 






(m) 




{Muj/ny/^ 2{Mojihyi^ 


■ (16) 



2MVtv 




QM^R 




(13) 


where S = U/U is the energy density and E is 
the total energy. In an external harmonic trap 
U(r) = iMw^jrp, the free energy can be written as 
F{T/Tp, kpV, kpR, hui/Ep) near resonance, and we find 

= (14) 

^ ' AMv ^ ^ AMR ^ 


where (V) denote the total potential energy due to the 
harmonic confinement. 

Momentum Distribution. The correlations encapsu¬ 
lated by the two contacts determine the tail of mo¬ 
mentum distribution, which can be measured using 
time-of-flight imaging [Tj. Theoretically, the momen¬ 
tum distribution can be obtained by Fourier trans¬ 
forming the single-particle density matrix / 9 i(r,r') = 
Af-i/d3r"(^t(r)^t(i.")^(i.")^(r')). While the leading 
term of momentum distribution is entirely determined by 
the internal structure of the two-body density matrix in 
the asymptotic regime [cf. Eq. 0 ], the sub-leading term 
depends on the distribution of center of mass momen¬ 
tum P. Assuming that the distribution of P is isotropic, 
namely, (P^) = (P^)/3 etc., we find that, in fact for 

kp k 1 /ro |TT] 
nic = E 

m 

StT^ dQOjq) P 

~ Z Up 

ATr^J2m.,p,3rip.j,m{JdqajqyP^ 

”^3 Ufc4 

which shows that while determines the 

strength of the leading I/fc^ [35j, the sub-leading term 
1 /P is not determined solely by the contact 
In fact, the angular dependences of the momentum dis¬ 
tribution is also not purely p-wave, but has an s-wave 
component. We also note that a sub-leading term in 
momentum distribution relating to the s-wave effective 
range is found in Ref. [10] . 

Radio-frequency Spectroscopy. The rf coupling FI^i = 
hklj-f f dr-ipl(r)ip(r) transfers fermions into an initially 
empty spin state |e), where Urf is the rf Rabi frequency. 
For a perturbative Urf, the transfer rate can be written 
as Frf(w) = (27T/h)J^. j pil(flH,fli}l^S(huj + Ei - Ef), 




MY 


327r2C'6 

MY 


|Ulm(fc)|" (15) 


Static Structure Factor. By definition, the static struc¬ 
ture factor S'(q) = 27r ^ ^ pi| (/| / Prn(r) exp(—iq • 
r)|i)p and can be measured by Bragg spectroscopy |iK] . 
Here n(r) is the density operator and other notations are 
the same as before. S'(q) can be obtained directly by 


Fourier transforming p 2 i Eq. (11), and diverges linearly 
in the limit q ^ oo. It is cut off by the short-range 
potential U(r) and will be limited by l/rg. 

Photo-association Spectroscopy. Photo-association 
has been used to measure the fraction of closed chan¬ 
nel molecules in two-component Fermi gases [46] . 
which is related to the s-wave contact [3 Ej- In 
the case of a p-wave resonance, if the internal wave 
function Pm(i’) of the relevant excited molecule has 
a specific projection m along the z direction, namely 
Pm(i’) Yimir), the transition rate is given by rp’r^(a;) = 

2nhnl,j:^jP,\{f\Jd^rd^Rgy{r)^liR)yM + 

r/2)i/;(R — r/2)|i)p<5(/ia; + Ei — Ef), with Upa the 
Rabi frequency and (/)l,j(R) the molecule creation 
operator. 

Since usually the final molecular state has a finite de¬ 
cay rate 7 , 6{Huj + Ei — Ef) in the expression of rp’r'^(w) 
should be replaced by a Lorentzian {hj/2)/[{huj + Ei — 
Ef)"^ A- {h'j/2)y. Typically 7 ^ 10 MHz much larger 
than the energy scales associated with the spatial mo¬ 
tion of the Fermi gas. As a result, when the to of the 
photo-association laser is tuned to resonance, h'y domi¬ 
nates over typical values of hui + Ei — Ef, and we can 
approximate the Lorentzian by 2/h'j, 


p(m) _ ^ff{m)Q2 

pa ^ ^ ^pa 


d^rgy{r)Yimif) 




(17) 


The Franck-Condon factor can be computed once Pm(r) 
is known. What is important here is that it depends 
only on two-body physics, so the many-body dependence 
is encapsulated in The contribution from 

is smaller by a factor (fcpro)^ if the excited molecular 
state is of extension cq. In the case when the photo¬ 
association process does not distinguish between hnal 
molecular states of different m, the total transition rate 
will be the sum of the individual Fp’T^ 

Virial Expansion for p-wave Contacts. At high tem¬ 
peratures, the effects of the interaction can be taken into 
account by the second virial expansion [47] . The change 
of the free energy of the spinless fermions 6F = F — Fq is 
given by SF/kpT = —2y/2NnX^b2, where Fq is the free 
energy without interactions and A = h/MjfMkpfT■ The 
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second virial coefficient is given by 
^2 = 3 


dk d5{k) ^ 

TT dk 


■ (18) 


Let Cv^R = by Ibe adiabatic theorems, 

§ = 8^/2™A^. ^ = 8^/2™A (19) 

When = 0, dh 2 ldv-^ = {3/V^)XR^K{X/{RV^)) 
with 


hviv) ='n + ri'^ 


dx 


(1 — e ^ )(77^ + 3x^) 


r(772a; +; 


,3'|2 


( 20 ) 


a.nddb 2 /dR ^ = (3/-\/27r)A/iij(A/(i?v^)) with 


hniv) = Erfc(77). 

a/tt 


( 21 ) 



FIG. 1. (Color online.) Dependence of the contacts Cv (the 
solid line) and Cr (the dashed line) on Eh/Ep for T/Tp = 2 
and their temperature dependences (inset) at a p-wave reso¬ 
nance from the virial expansion. We have taken kpR = 1/25. 


Figure shows the dependences of Cy and Cr as a func¬ 
tion of Eb/Ep for T/Tp = 2 and kpR = 1/25, appropriate 
for the case of with a typical density of 2 x 10^® m“^ 
at the p-wave resonances near B = 198.5G HU. Here we 
note that while Cy decreases monotonically as —Eb/Ep 
increases, Cr shows non-monotonic behavior and reaches 
a maximum when —Eb/Ep ~ 2, where it is comparable to 
Cy, if non-dimensionalized by kp (see Fig. [^. The tem¬ 
perature dependence of the contacts Cy and Cr at v~^ = 
0 is shown in the inset, for which Cr is much smaller 
than Cy] the magnitude of both grows with increasing 
R. In the ternperature regime Tp < T h^/2MR^, 
Eqs. @ and @ give Cy « 12\/2iVnA3i? ~ T"®/® and 
Cr « 247r7Vn.R^l — iR^MkpT which should be 
contrasted with T“^-dependence for s-wave contact [45] . 

Away from p-wave resonances where the scattering vol¬ 
ume V is small, Eq. (18) gives Cy = 12tt'^ Nnv"^ / X^ and 
Cr = 3607r®A^nu^/A^ when contribution from the deeply 
bound state is excluded. In this limit, the scaling Cy ~ 
and Cr ^ is also expected from perturbation calcula¬ 
tions for small v |54|, which indicates the irrelevance of 
Cy and Cr, in Fermi gases close to an s-wave resonance. 

Discussion. Our derivation of the p-wave contacts is 
based on the single-channel model which does not take 
into account explicitly the presence of closed-channel 
molecules, as in the case of a Feshbach resonance. The 
same results shall be obtained for a two-channel model, 
provided that the closed-channel molecule is small (com¬ 
parable to ro), which is typically the case. This is because 
all the arguments so far depend only on the properties of 
the two-body wave function or two-body density matrix 
in the asymptotic regime, which, in our derivation, de¬ 
pend only on the scattering volume v and effective range 
R, irrespective of whether they arise from a shape reso¬ 
nance or a Feshbach resonance. For actual atomic sys¬ 
tems, the van der Waals potential modifies the p-wave 
scattering phase shift by introducing a term a/k^ in the 


effective range expansion [5^ |4U. However, close to a 
Feshbach resonance, it was shown that a ~ 1/u^, whose 
effects are thus negligible |26l[50]. As a result, we expect 
that our main results Eqs. (Ill to 0 to remain true 
close to a p-wave Feshbach resonance. 

Resonances for different \m\ can be split due to mag¬ 
netic dipole-dipole couplings [501. For the m = 
0 and m = ±I resonances around B = 198.5 G 
are split by about 0.5 G [20|- To take this into ac¬ 
count, we can introduce phase shifts for different m, 
coti5^™^ = —1/v^'^^k^ — 1/R^^^k. Likewise, we can es¬ 
tablish the relation dE/d{l/v^'^'^) = —h'^Cy^'^/2M and 
d£;/d(l/R(™)) = I2M while Eqs. ^ and ^ 

to ( [l7| stay intact. 
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Note added. During the final preparation of this 
manuscript, closely related work by Yoshida and Ueda 
appeared |5T], in which they discuss one of the contacts, 
Cy, using a two-channel model. 


[1] W. Zwerger, ed., The BCS-BEC Crossover and the Uni¬ 
tary Fermi Gas, (Springer-Verlag 2011). 

[2] Tin-Lun Ho, Phys. Rev. Lett. 92 090402 (2004). 

[3] Shina Tan, Ann. Phys. N.Y. 323, 2952 (2008). 
























5 


[4] E. Braaten and L. Platter, Phys. Rev. Lett. 100, 205301 
(2008). 

[5] S. Zhang and A. J. Leggett, Phys. Rev. A 79, 023601 
(2009). 

[6] F. Werner, L. Tarruell, Y. Castin, Enro. Phys. J. B 68, 
410 (2009). 

[7] J. T. Stewart, J. P. Gaebler, T. E. Drake, and D. S. Jin, 
Phys. Rev. Lett. 104, 235301 (2010). 

[8] Y. Sagi, T. E. Drake, R. Paudel, and D. S. Jin, Phys. 
Rev. Lett. 109, 220402, (2012). 

[9] S. Hoinka, M. Lingham, K. Fenech, H. Hn, C. J. Vale, J. 
E. Drut, and S. Gandolfi, Phys. Rev. Lett. 110, 055305 
(2013). 

[10] M. Valiente, N.T. Zinner and K. M0lmer, Phys. Rev. 84, 
063626 (2011). 

[11] M. Valiente, N.T. Zinner and K. M0lmer, Phys. Rev. A 
86, 043616 (2012). 

[12] G. A. Regal, C. Ticknor, J. L. Bohn, and D. S. Jin, Phys. 
Rev. Lett. 90, 053201 (2003). 

[13] J. Zhang, E. G. M. van Kempen, T. Bourdel, L. 
Khaykovich, J. Cubizolles, F. Chevy, M. Teichmann, L. 
Tarruell, S. J. J. M. F. Kokkelmans, and C. Salomon, 
Phys. Rev. A 70, 030702(R) (2004). 

[14] Kenneth Gunter, Thilo Stoferle, Henning Moritz, 
Michael Kohl, and Tilman Esslinger, Phys. Rev. Lett. 
95, 230401 (2005). 

[15] C. H. Schunck, M. W. Zwierlein, C. A. Stan, S. M. F. 
Raupach, W. Ketterle, A. Simoni, E. Tiesinga, C. J. 
Williams, and P. S. Julienne, Phys. Rev. A 71, 045601 
(2005). 

[16] J. P. Gaebler, J. T. Stewart, J. L. Bohn, and D. S. Jin, 
Phys. Rev. Lett. 98, 200403 (2007). 

[17] J. Fuchs, C. Ticknor, P. Dyke, G. Veeravalli, E. Kuhnle, 
W. Rowlands, P. Hannaford, and C. J. Vale, Phys. Rev. 
A 77, 053616 (2008). 

[18] Y. Inada, M. Horikoshi, S. Nakajima, M. Kuwata- 
Gonokami, M. Ueda, and T. Mukaiyama, Phys. Rev. 
Lett. 101, 100401 (2008). 

[19] Takuya Nakasuji, Jun Yoshida, and Takashi Mukaiyama, 
Phys. Rev. A 88 012710 (2013). 

[20] C. Ticknor, C. A. Regal, D. S. Jin, and J. L. Bohn, Phys. 
Rev. A 69, 042712 (2004). 

[21] F. Chevy, E.G.M. van Kempen, T. Bourdel, J. Zhang, 
L. Khaykovich, M. Teichmann, L. Tarruell, S.J.J.M.F. 
Kokkelmans, C. Salomon, Phys. Rev. A 71, 062710 
(2005). 

[22] J. Levinsen, N. R. Cooper, and V. Gurarie, Phys. Rev. 
Lett. 99, 210402 (2007). 

[23] J. Levinsen, N. R. Cooper, and V. Gurarie, Phys. Rev. 
A 78, 063616 (2008). 

[24] L. Pricoupenko, Phys. Rev. Lett. 96, 050401 (2006). 


[25] M. Jona-Lasinio, L. Pricoupenko, and Y. Castin, Phys. 
Rev. A 77, 043611 (2008). 

[26] P. Zhang, P. Naidon, and M. Ueda, Phys. Rev. A 82, 
062712 (2010). 

[27] E. Braaten, P. Hagen, H.-W. Hammer, and L. Platter, 
Phys. Rev. A 86, 012711 (2012). 

[28] Y. Nishida, S. Moroz, and D. T. Son, Phys. Rev. Lett. 
110, 235301 (2013). 

[29] A. G. Volosniev, D. V. Fedorov, A. S. Jensen, N. T. Zin¬ 
ner, J. Phys. B 47, 185302 (2014). 

[30] S. Moroz and Y. Nishida, Phys. Rev. A 90, 063631 
(2014). 

[31] C. Gao, J. Wang and Z. Yu, arXiv:1412.3566, 

[32] S. G. Peng, S. Tan and K. Jiang, Phys. Rev. Lett. 112, 
250401 (2014). 

[33] D.K. Gridnev, J. Phys. A: Math. Theor. 47 505204 
(2014) 

[34] T. Y. Gao, S. G. Peng, and K. Jiang, Phys. Rev. A 91, 
043622 (2015). 

[35] Y. Ohashi, Phys. Rev. Lett. 94, 050403 (2005). 

[36] V. Gurarie, L. Radzihovsky, Annals of Physics 322, 2 
(2007). 

[37] K. B. Gubbels and H. T. C. Stoof, Phys. Rev. Lett. 99, 
190406 (2007). 

[38] D. Inotani, R. Watanabe, M. Sigrist, and Y. Ohashi, 
Phys. Rev. A 85, 053628 (2012). 

[39] C. Luciuk, S. Trotzky, S. Smale, Zhenhua Yu, Shizhong 
Zhang and J. H. Thywissen, arXiv: 1505.08151 (2015). 

[40] F. Werner and Y. Castin, Phys. Rev. A, 86, 013626 

( 2012 ). 

[41] See supplementary materials for derivations of the shal¬ 
low p-wave two-body bound states, derivation of equa¬ 
tions (3) and (4), and the tail of momentum distribution. 

[42] M. Punk and W. Zwerger, Phys. Rev. Lett. 99, 170404 
(2007). 

[43] W. Schneider and M. Randeria, Phys. Rev. A 81, 021601 

( 2010 ). 

[44] E. Braaten, D. Kang, L. Platter, Phys. Rev. Lett. 104, 
223004 (2010) 

[45] S. Hoinka, M. Lingham, K. Fenech, H. Hu, C.J. Vale, 
J.E. Drut, and S. Gandolfi, Phys. Rev. Lett. 110, 055305 
(2013). 

[46] G. B. Partridge, K. E. Strecker, R. I. Kamar, M. W. Jack, 
and R. G. Hulet, Phys. Rev. Lett. 95, 020404 (2005). 

[47] T-L Ho and E.J. Mueller, Phys. Rev. Lett. 92, 160404 
(2004). 

]48] Z. Yu, G.M. Bruun, and G. Baym, Phys. Rev. A 80, 
023615 (2009). 

[49] B. Gao, Phys. Rev. A 58, 4222 (1998). 

[50] B. Gao, Phys. Rev. A 84, 022706 (2011). 

[51] Shuhei M. Yoshida and Masahito Ueda, Phys. Rev. Lett. 
115, 135303 (2015) 


